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What happened to the central cores of tidally destructed dark matter clumps in the Galactic halo? 
We calculate the probability of surviving of the remnants of dark matter clumps in the Galaxy by 
modelling the tidal destruction of the small-scale clumps. It is demonstrated that a substantial 
fraction of clump remnants may survive through the tidal destruction during the lifetime of the 
Galaxy if the radius of a core is rather small. The resulting mass spectrum of survived clumps is 
extended down to the mass of the core of the cosmologically produced clumps with a minimal mass. 
Since the annihilation signal is dominated by the dense part of the core, destruction of the outer 
part of the clump affects the annihilation rate relatively weakly and the survived dense remnants 
of tidally destructed clumps provide a large contribution to the annihilation signal in the Galaxy. 
The uncertainties in minimal clump mass resulting from the uncertainties in neutralino models are 
discussed. 

PACS numbers: 12.60.Jv, 95.35.+d, 98.35.Gi 



I. INTRODUCTION 

According to current observations, about 30% of the 
mass of the Universe is in a form of cold dark matter 
(DM). The nature of DM particles is still unknown. The 
cold DM component is gravitationally unstable and is ex- 
pected to form the gravitationally bounded clumpy struc- 
tures from the scale of the superclusters of galaxies and 
down to very small clumps of DM. The large-scale DM 
structures are observed as the galactic halos and clus- 
ters of galaxies. They are also seen in numerical simu- 
lations. Theoretical study of DM clumps are important 
for understanding the properties of DM particles because 
annihilation of DM particles in small dense clumps may 
result in visible signal. The DM clumps in the Galaxy 
can produce the bright spots in the sky in the gamma or 
X-bands [l[ . A local annihilation rate is proportional to 
the square of the DM particle number density. Thus, the 
annihilation signal from small clumps can dominate over 
diffuse component of DM in the halo. 

The cosmological formation and evolution of small- 
scale DM clumps have been studied in numerous works 
[1, S, II H B 0,11 S EH EH, El, El • The minimum mass 
of clumps (the cutoff of the mass spectrum), M min is de- 
termined by the collisional and collisionlcss damping pro- 
cesses (see e. g. @ and references therein). Recent calcu- 
lations show that the cutoff mass is related to the fric- 
tion between DM particles and cosmic plasma similar to 
the Silk damping. In the case of the Harrison- Zeldovich 
spectrum of primordial fluctuations with CMB normal- 
ization, the first small-scale DM clumps are formed at 
redshift z ~ 60 (for 2a fluctuations) with a mean density 
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7 x 10~ 22 g cm -3 , virial radius 6 x 10 -3 pc and internal 
velocity dispersion 80 cm s _1 respectively. Only a very 
small fraction of these clumps survives the early stage 
of tidal destruction during the hierarchial clustering [4|] . 
Nevertheless, these survived clumps may provide the ma- 
jor contribution to the annihilation signal in the Galaxy 
0, 0, EH EE EH- At a high redshift, neutralinos, con- 
sidered as DM particles, may cause the efficient heating 
of the diffuse gas [ItJ due to annihilation in the dense 
clumps. 

One of the unresolved problem of DM clumps is a value 
of the central density or core radius. Numerical simula- 
tions give a nearly power density profile of DM clumps. 
Both the Navarro-Frenk- White and Moore profiles give 
formally a divergent density in the clump center. A the- 
oretical modelling of the clump formation [l8j predicts a 
power-law profile of the internal density of clumps 

**(!■) = — P {-) , (1) 

where p and R are the mean internal density and a radius 
of clump, respectively, f3 ~ 1.8 — 2 and Pi n t( r ) = at r > 
R. A near isothermal power-law profile ([T]) with (3 ~ 2 
has been recently obtained in numerical simulations of 
small-scale clump formation [19j |. 

It must be noted that density profiles of small-scale 
DM clumps and large-scale DM haloes may be differ- 
ent. The galactic halos are well approximated by the 
Navarro-Frenk- White profile outside of the central core 
where dynamical resolution of numerical simulations be- 
comes insufficient. Different physical mechanisms are en- 
gaged for formation of a central core during the formation 
and evolution of clumps. A theoretical estimation of the 
relative core radius of a DM clump x c — R c /R was ob- 
tained in [l8[ from energy criterion, x c = R c /R — 8l q , 
where 8 eq is a value of density fluctuation at the begin- 
ning of a matter-dominated stage. A similar estimate 
for DM clumps with the minimal mass ~ 10 -6 M Q orig- 
inated from 2a fluctuation peaks gives S eq ~ 0.013 and 
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R c /R ~ 1.8 x 10~ 5 , respectively. In the core radius 
x c ~ 0.3f -2 has been obtained, where v is a relative 
height of the fluctuation density peak in units of dis- 
persion at the time of energy-matter equality (see also 
Section [V}. This value is a result of the influence of tidal 
forces on the motion of DM particles in the clump at the 
stage of formation. This estimate may be considered as 
an upper limit for the core radius or as the break-scale 
in the density profile, e.g., a characteristic scale in the 
Navarro-Frenk- White profile. It could be that a real core 
radius, where the density ceases to grow, is determined 
by the relaxation of small-scale perturbations inside the 
forming clump [20]. Another mechanism for core forma- 
tion arises in the "meta-cold dark matter model" due to 
late decay of cold thermal relics into lighter nonrelativis- 
tic particles with low phase-space density [2ll. |22|. 

Nowadays, numerical simulations have a rather low 
space resolution in the central region of clumps to de- 
termine the core radius. The only example with some 
indication to presence of a core with radius x c ~ 10~ 2 is 
numerical simulation of small-scale clump formation . 
The special numerical simulations with a sufficiently high 
space resolution to reveal the real core radius are very re- 
quested. 

In this work, we consider the relative core radius 
x c = R c /R of DM clumps as a free parameter in the 
range 0.001 — 0.1. We investigate the dependence of the 
probability of clump survival in the Galaxy on this pa- 
rameter under the action of tidal forces from galactic disk 
and stars. As a preferred, value we consider x c ~ 10~ 2 in 
spite of the numerical simulations (l9j . The correspond- 
ing annihilation rate of DM particles is proportional to 
their squared number density, and thus is very sensitive 
to the value of the core radius. 

In our earlier works @, fl2| , we used a simplified cri- 
terium for a tidal destruction of clump. Namely, we pos- 
tulated that the clump is destructed if a total tidal energy 
gain after several disk crossings (or collisions 

with stars) becomes of order of initial binding energy of 
a clump \E\, i.e. 



\E\ 



(2) 



where summation goes over the successive disk crossings 
(or encounters with stars). This criterium is justified 
in the cosmological context of the DM clump formation 
because both the formation of the density profile of the 
clump and its tidal heating proceed during the same time 
of nonlinear evolution of density perturbation. For the 
Galaxy case, a more detailed consideration is needed to 
describe a tidal destruction of DM clumps by stars. An 
improved ap proach includes a gradual mass loss of sys- 
te ms [23l [2J, |25| , in particular, by small-scale DM clumps 

I, Hi 

In this work, we will describe a gradual mass loss of 
small-scale DM clumps assuming that only the outer lay- 
ers of clumps are involved and influenced by the tidal 
stripping. Additionally, we assume that inner layers of a 



clump are not affected by tidal forces. In this approx- 
imation, we calculate a continuous diminishing of the 
clump mass and radius during the successive Galactic 
disk crossings and encounters with the stars. We accept 
now for criterium of clump destruction the diminishing 
of the radius of tidally stripped clump down to the core 
radius. An effective time of mass loss for the DM clump 
remains nearly the same as in our previous calculations 
|12j. However, the clump destruction time has now quite 
different physical meaning: it provides now a character- 
istic time-scale for the diminishing of the clump mass 
and size instead of the total clump destruction. This 
means that small remnants of clumps may survive in the 
Galaxy. Respectively, these remnants would be an ad- 
ditional source of amplification of the DM annihilation 
signal in the Galaxy. 



II. TIDAL DESTRUCTION OF CLUMPS BY 
DISK 

The kinetic energy gain of a DM particle with respect 
to the center of a clump after one crossing of the Galactic 
disk is 
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where m is a constituent DM particle mass, Az is a ver- 
tical distance (orthogonal to the disk plane) of a DM 
particle with respect to the center of the clump, i> z . c is 
a vertical velocity of the clump with respect to the disk 
plane at the moment of disk crossing, and A(a) is the 
adiabatic correction factor. A gravitational acceleration 
near the disk plane is 



g m (r) = 2nGa s (r), 



(4) 



where we use an exponential model for a surface density 
of disk 



as{r) = ^L e -r/r 



2^r 2 



(5) 



with M d = 8x 10 10 M Q , r = 4.5 kpc. 

The factor A{a) in §5§ describes the adiabatic protec- 
tion from slow tidal effects [28j. This adiabatic correc- 
tion, further on referred to as the Weinberg correction, is 
defined as an additional factor A(a) to the values of en- 
ergy gain in the momentum approximation. This factor 
satisfies the following asymptotic conditions: A(a) = 1 
for a < 1 and A(a) -C 1 for a > 1. In [23|. the following 
fitting formula was proposed: 



A(a) = (l + a 2 )- 3/2 . 



(6) 



Here the adiabatic parameter a — lut^, where u> is an or- 
bital frequency of DM particle in the clump, ~ Hd/v ZtC 
is an effective duration of gravitational tidal shock pro- 
duced by the disk with a half-thickness Ed- For tidal 
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interactions of clumps with stars in the bulge and the 
halo the duration of the gravitational shock can be es- 
timated as t s ~ l/v Ie i, where I is an impact parameter 
and v re i is a relative velocity of a clump with respect to 
a star. 

As a representative example, we consider the isother- 
mal internal density profile of a DM clump 



PintO) 
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with a cutoff at the virial radius R: p(r) — at 
r > R. A corresponding mass profile of a clump is 
M(r) = Mi(r/R), where Mi is an initial mass of a clump 
at the epoch of the Galaxy formation. With this mass 
distribution, a circular velocity inside a clump is indepen- 
dent of the radius, v rot = (GM(r) /r) 1 / 2 = {GM./R) 1 ' 2 . 
A gravitational potential corresponding to the density 
profile ([7]) is 0(r) = v^ ot [log (r/R) — 1]. Let us define a di- 
mensionless energy of the DM particle e = E/(mv 2 ot ) and 
gravitational potential tp(r) — <f)(r)/v 2 ot — h\{r/R) — 1. 
An internal density profile Pxotir) and the distribution 
function of DM particles in the clump f c i(s) are related 
by the integral relation (29j 
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The corresponding isothermal distribution function is 
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Note what this distribution function provides only an 
approximate representation of ([7]), far from the cutoff 
radius R. Nevertheless, this approximation is enough for 
our estimates of tidal destruction of the DM clumps. 

By using the hypothesis of a tidal stripping of outer 
layers of a DM clump, we see that a tidal energy gain 
5s causes the stripping of particles with energies in the 
range — Se < e < 0. A corresponding variation of density 
at radius r is 



Sp(r) = 2 5 / 2 TT / ^e-yj(r)f cl (e)de. 



(10) 
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In this equation, the tidal energy gain (J3J by different DM 



particles is averaged over angles, so as ((A* 
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A resulting total mass loss by a DM clump during one 
crossing of the Galactic disk is 



S M = -Am I r 5p{r) dr. 
Jo 

Let us specify the dimensionless quantities 

Qd-7, 2 r <- i b d =—GpiT d , 



(11) 



(12) 



where pi = 3Mi/(47ri? 3 ) is a initial mean density of 
clump. For the most parts of clumps Qd <C 1 with a 
typical value Qd ~ 0.03. In the limiting case Qd <C 1 
and in the absence of the adiabatic correction, Sd = 0, 
the integrals (fT0|) can be calculated analytically. In a 
general case, the fitting formula for the mass loss of a 
clump during one passage through the Galactic disk is 



(SM\ 



-0.13Q d exp (-1.58^ 
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(13) 



Now, we calculate the tidal mass loss by clumps using 
a realistic distribution of their orbits in the the halo. 
The method of calculation is similar to the one used in 
[l2l |. but instead of the rough energetic criterium for a 
tidal clump destruction ^ we will assume now a gradual 
decreasing of clump mass and size. 

Let us choose some particular clump moving in the 
spherical halo with an orbital "inclination" angle 7 be- 
tween the normal vectors of the disk plane and orbit 
plane. The orbit angular velocity at a distance r from 
the Galactic center is d<fi/dt = J/(mr 2 ), where J is an 
orbital angular momentum of a clump. A vertical veloc- 
ity of a clump crossing the disk is 



J 

mr s 



sin 7, 



(14) 



where r s is a radial distance of a crossing point from 
the Galaxy center. There are two crossing points (with 
different values of r s ) during an orbital period. 

The standard Navarro-Frenk- White profile of the DM 
Galactic halo is 



Po 



(r/L)(l + r/LY 



where L — 45 kpc, po — 5 x 1O 6 M kpc 3 . 
introduce the dimensionless variables: 



x = — , pn(x) 
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(15) 

It useful to 
(16) 
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where $0 = —^nGpoL 2 , E or \> is a total orbital energy of 
a clump. With these variables, the density profile of the 
halo (fT5")) is written as 



Ph(x) = 



1 



x(l + x) 2 ' 



(18) 



A gravitational potential ipix), corresponding to density 
profile (HU) is 



tp(x) 



1 



log(l + x) 



(19) 



An equation for orbital turning points, r 2 — 0, for DM 
clumps in the potential (fH?)) is 



1 



log(l + x) 



y 



(20) 
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From pop , one can find numerically the minimum x m i n 
and maximum x max radial distance of a clump from the 
Galactic center as a function of orbital energy e and 
square of angular momentum y. Denoting p — cos 9, 
where 9 is an angle between the radius- vector f and the 
orbital velocity v, we have y = (1 — p 2 )x 2 [e — ip(x)). As 
we assumed above, the unit vectors v/v are distributed 
isotropically at each point x, and, therefore, p has a uni- 
form distribution in the interval [—1, 1]. 

The relation between the density profile Ph{ x ) and the 
distribution function is given by the same equation ([8]) 
with an obvious substitution / c i => F(e), where the dis- 
tribution function F(e) for a halo profile (fT8|) can be fit- 
ted as m 



3/2-5/2 



ln(l - e) 



Here F 1 = 9.1968 x 10" 



(21) 



-2.7419, P = EP*(l-e) i . 



(Pl,P2,P3,Pi) 

An inte 
back is 

T c (x,e,p) 



(0.3620, -0.5639, -0.0859, -0.4912). 



ds 



V27rGpo 7 - ^(s) - y/s 2 
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An angle of orbital precession during the time T c /2 is 

^max 

ds 



4> = v 
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V e V'(s) - y/s 2 



- 7T < 0. 



(23) 



Therefore, an orbital period is longer than T c and is given 
by 

T t = T c (l + 0/tt) \ (24) 

Choosing a time interval AT much longer than a clump 
orbital period T t , but much shorter than the age of the 
Galaxy to, i.e., T t <C AT <C to, we define an averaged 
rate of mass loss by a selected clump under influence of 
tidal shocks in successive disk crossings 
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' dM\ 

~dT) d 

where (SM/M)d is given by (fill)) and summation goes 
over all successive crossing points (odd and even) of a 
clump orbit with the Galactic disk during the time inter- 
val AT. According to (dJ) and (TT4"]) the g m and u ZiC both 
depend on the radius x — r/L. One simplification in cal- 
culation of ([25]) follows from the fact that a velocity of 
orbit precession is constant. For this reason the points of 
successive odd crossings are separated by the same an- 
gles (j) from (|2"3"j) . The same is also true for successive 
even crossings. Using this simplification, we transform 
the summation in (|25p to integration: 



1 ^ r sM 

AT^ l~M" 
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, dx 



where 



,1/2 



dx 



(26) 



is an equation for the clump orbit in the halo. The 
method described will be used in the Section IIVI for the 
final calculations. 



III. TIDAL DESTRUCTION OF CLUMPS BY 
STARS 

Now, we calculate the diminishing of a clump mass 
due to a tidal heating by stars in the Galaxy by using 
the same hypothesis of the the preferable stripping of the 
outer clump layers. During a single close encounter of a 
DM clump with a star, the energy gain of a constituent 
DM particle in the clump with respect to clump center 

is a 
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where to* is a star mass, I is an impact parameter, v ro \ 
is a relative star velocity with respect to a clump, Az = 
rcosijj, r is a radial distance of a DM particle from the 
clump center, and ip is an angle between the directions 
from the clump center to the DM particle and to the 
point of the closest approach of a star. Using the same 
method as in Sec. HH we calculate a relative mass loss 
by the clump (SM/M) S during a single encounter with a 
star and obtain the same fitting formula as (|13[) but with 
substituting the dimensionless parameters, Qd Q s and 
Sd => S s , where 
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(28) 



where t s ~ l/v Ie \. 

A DM clump acquires the maximum energy gain dur- 
ing a single encounter with a star when the impact pa- 
rameter I ~ R. Using the relation 



dt 



dx 



2^/2irGpo y/s — ip(x) — y/x 2 ' 



(29) 



and integrating over all impact parameters I > R, we 
calculate an averaged rate of mass loss by a clump during 
successive encounters with stars 



1 fdM 
M \~~dT 



(30) 



oo x„ 



2T tv /2 



1 f 7 dsnMvrei fSM\ 

2-kGpq J J ^e-V(s)-2//s 2 V M A 



R %n 



where n*(r) is a radial number density distribution of 
stars in the bulge and halo. A DM clump moves through 
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the medium with a varying value of along the clump 
orbit. In contrast to the case of the disk crossing, the 
precession of the clump orbit during an orbital period 
does not influence the mass loss due to encounters with 
stars. Additionally, the mass loss due to encounters with 
stars is independent of the inclination of clump orbits in 
the case of a spherically symmetric distribution of stars 
in the bulge and halo. 

Using the results of [3l| . we approximate the radial 
number density distribution of stars in the bulge in the 
radial range r = 1 — 3 kpc as 



nb,*{r) = (p h /m*)exp \-{r/n) 



1.61 



(31) 



where pb = 8M /pc 3 and n, = 1 kpc. A corresponding 
number density distribution of the halo stars at r > 3 kpc 
outside the Galactic plane can be approximated as 



n h ,*{r) = ( / Oh/TO*)(ro/r) 3 



(32) 



where = 0.4M Q and r© = 8.5 kpc. According to 
[32| , in the region between r = 1 and 40 kpc a total mass 
of stars is 4 x 10 8 Af Q with a star density profile oc r~ 3 . 
These data correspond to ph — 1.4 x 10~ 5 M /pc 3 in 
(|3"2")l . We neglect in our calculations the oblationes of the 
stellar halo 13211. 



IV. SURVIVING FRACTION OF CLUMPS 

From Eq. (|3|) it is seen that the tidal forces influence 
mainly the outer part of the clump (where Az is rather 
large). Further, we will use our basic assumption that 
only outer layers of a clump undergo the tidal stripping, 
while the inner parts of a clump are unaffected by tidal 
forces. Thus, we assume that a clump mass M = M(t) 
and radius R = R(t) are both diminishing in time due to 
the tidal stripping of outer layers, but its internal density 
profile remains the same as given by Eq. ([7]), e.g., for 
the isothermal density profile M(t) oc R(t) and p(t) oc 
M(t)~ 2 . Combining together the rates of mass loss ([23]) 
and (j30|) due to the tidal stripping of a clump by the disk 
and stars respectively, we obtain the evolution equation 
for a clump mass: 



dM _ fdM\ fdM 
~~dT ~ \ dT)d + 



(33) 



In the following, we solve this equation numerically start- 
ing from the time of Galaxy formation at to — t(j up to 
the present moment to- In numerical calculations, it is 
convenient to use the dimensionless variables: t/to for 
time and M/Mi for a clump mass, where Mi is an initial 
clump mass. The adiabatic correction provides generally 
only a small effect. In the absence of adiabatic correction 
or, equivalently, at Sd — S s ~ the evolution equation 
(|3"3")l has a simple form 



dp 
~dt 



P_ 

t s 



p_ 

t d 



(34) 



where p = M(t)/Mi and parameters td and t s are inde- 
pendent of p. The solution of this equation 



2td 

{2t d + t s )exp(2t /t s )-t, 



(35) 



represents a good approximation to numerical solution of 
(f3"3")) with the adiabatic correction taken into account. 

The most important astrophysical manifestation of 
DM clumps is a possible annihilation of constituent DM 
particles. The crucial point is a dominance of the cen- 
tral core of a clump in annihilation signal if clumps have 
a steep enough density profile. Namely, annihilation of 
DM particles in a clump core will prevail in a total anni- 
hilation rate in a single clump with a power-law density 
profile |U) if /3 > 3/2 and x c = R c /R < 1. More specif- 
ically, the quantity N oc £ 4:Trr' 2 dr' pf nt (r') practically 
does not depend on r, if r 3> ro- As a result, the anni- 
hilation luminosity of a DM clump with approximately 
an isothermal density profile (/3 ~ 2) will be nearly con- 
stant under influence of tidal stripping until a clump ra- 
dius diminishes to its core radius. In other words, in 
the nowadays Galaxy the remnants of tidally stripped 
clumps with x c < p(to) <C I, where p(t) — M(t)/Mi 
and to ~ 10 10 yrs is the Galaxy age, obeys the evolution 
equation (|33p and have the same annihilation luminosity 
as their progenitors with p = 1. 

By using the evolution equation (|33[) . we now calculate 
the probability P of the survival of clump remnant during 
the lifetime of the Galaxy. Let us choose some arbitrary 
point in the halo with a radius- vector r and an angle a 
with a polar axis of the Galactic disk. Only the clump 
orbits with an inclination angle 7r/2 — a < "f < ir /2 pass 
through this point. A survival probability for clumps can 
be written now in the following form 



1 sin a 
p, s 47TV2 f.f, 

P{x,a) = —— dp a cos 7 

p(x)sma 



(36) 







x / de[e-^{x)}^ 2 F{e)e[p(to)~x c }. 

ip(x) 

In this equation, p(x) is a density profile of the halo from 
(HHJ), p = cos 8, 9 is an angle between the radius- vector 
r and the orbital velocity of clump, is the Heaviside 
function, ij)(x) is the halo gravitational potential from 
(fT5|) . F(e) is a distribution function of clumps in the halo 
from (|2"Tj) , p(to) depends on all variables of the integra- 
tion, and x c = R c /R is an initial value of the clump core. 
The function p(to) is calculated from numerical solution 
of evolution equation (|33[) . If p(to) > x c , the clump rem- 
nant is survived through the tidal destruction by both 
the disk and stars. The annihilation rate in this remnant 
would be the same as in the initial clump. On the con- 
trary, in the opposite case, when p(tc) < x c , the clump 
is totally destructed because (i) the core is not a dynam- 
ically separated system and composed of particles with 



6 




P(r) 



r, kpc 



FIG. 1: The survival probability P(r,p) plotted as a function 
of distance from the Galactic center r and a mean internal 
clump density p in the case x c — 0.1. It gives the normalized 
fraction of DM clumps in the halo P calculated from (|36|) . 
which survives the tidal destruction by the stellar disk and 
the halo stars. 




P(r) 



r, kpc 



FIG. 2: The same as Fig. [T] but for the case x c = 0.05. 



initial radius of the clump Ri depends only on the one pa- 
rameter — the initial clump density, which also depends 
only on the initial clump mass Mj. 

The crucial result of numerical calculation of a sur- 
vival probability (|36|) for clumps with x c <C 0.05 is that 
P(x, a) ~ 1 everywhere. Even inside the bulge there 
are clumps which flying through the bulge from external 
regions. These means that clump remnants are mostly 
survived through the tidal destruction in the Galaxy. A 
noticeable diminishing P(x, a) < 1 near the center of 
Galaxy becomes apparent for clumps with x c > 0.05. It 
is understandable because with x c — > 1, we return to the 
previous criterium of tidal destruction of clumps and 
to a corresponding results for survival probability p], [12] • 
The survival probability P(r, a) numerically calculated 
from (|36p for the cases x c ~ 0.05 and 0.1 is shown in 
the Figs. [1] and O The dependence on a (an angle be- 
tween a radius- vector r and a polar axis of the Galactic 
disk) is very weak as it was shown in [l2T |. For this rea- 
son, we present the results only for an intermediate value 
a = 7r/4. The density of clumps is normalized to the 
density 7.3 x 10~ 23 g cm -3 valid for clumps with mass 
AI = 10 _6 M© originated from 2a density peaks in the 
case of power-law index of primordial spectrum of per- 
turbations n p = 1. 

It is worth to note that a tidal radius of a clump in the 
bulge is (33] ] 



^7, 



GM{r t ) 
d 2 <f>/dP 



(37) 



where uj p ~ [GMb(l) /P] 1 / 2 is an angular velocity at the 
pericenter (we consider here a circular orbit for simplic- 
ity), Mb(l) is a mass profile of the bulge and <f>(l) is a 
gravitational potential of the bulge. For the considered 
small-scale clumps r t > 0.2i?j, and, therefore, the tidal 
radius is not a crucial factor for destruction of clumps. 



extended orbits, and because (ii) a nearly homogeneous 
core is destructed easier than a similar object with the 
same mass but with a near isothermal density profile. 

We consider the small-scale DM clumps in the initial 
mass interval Mj = [10 _6 M©, 1 M©] originated from the 
2(7 (i.e. v = 2) peaks in the Harrison- Zeldovich pertur- 
bation spectrum. A reason is as follows. The DM clumps 
originated from initial density perturbations with v < 1 
were almost completely destructed by tidal interactions 
during the early stage of hierarchical clustering as it can 
be seen from the distribution function of clumps (I44|) (see 
below). On the contrary, the most dense DM clumps 
with v > 3 are mostly survived the stage of hierarchical 
clustering, but their number according to (|44[) is expo- 
nentially falls with v and is small. For this reason, we 
will use the following approximation: the DM clumps 
were originated on average from v ~ 2 peaks, and for 
any given mass M we do not consider the distribution of 
clumps over their densities. In this approximation, the 



V. COSMOLOGICAL DISTRIBUTION 
FUNCTION OF CLUMPS 

In this section, we provide calculations of a mass func- 
tion for the small-scale clumps in the Galactic halo by 
more transparent method than in our previous works 

The first gravitationally bound objects in the Universe 
are the DM clumps of minimum mass M m j n . A numeri- 
cal value of M m i n depends strongly on the nature of DM 
particle. Even in the case of a particular DM particle, 
e.g., neutralino, the calculated value of M m ; n can differ 
by many orders of magnitude for different sets of param- 
eters in mSUGRA model, see the Appendix IA II The 
larger scale clumps are formed later. The larger scale 
clumps host the smaller ones and are hosted themselves 
by next larger clumps. Major parts of small-scale clumps 
are destroyed by the tidal gravitational fields of their host 
clumps. At small-mass scales, the hierarchial clustering is 



7 



a fast and complicated nonlinear process. The formation 
of new clumps and their capturing by the larger ones are 
nearly simultaneous processes because at small-scales an 
effective index of the density perturbation power spec- 
trum is very close to a critical value, n — ► —3. The 
DM clumps are not totally virialized when they are cap- 
tured by hosts. The adiabatic invariants cannot prevent 
the survival of cores at this stage because there are not 
enough time for the formation of the singular density pro- 
files in clumps (an internal dynamical time of clump is of 
the same order as its capture time by a host). We use a 
simplified model to take into account the most important 
features of hierarchial clustering. 

In the model of spherical collapse (see for example 
[34]]), a formation time for a clump with an internal den- 
sity p is t = (K,p cq / p) 1/2 t cq , where n = 187r 2 , p cq = 
po(l + z cq) 3 is a cosmological density at the time of 
matter-radiation equality t eq , 1 + z cq = 2.35 x 10 4 fi OT /i 2 , 
and po = 1.9 x 10~ 29 f2 m /i 2 g cm~ 3 . The index "eq" refers 
to quantities at the time of matter-radiation equality t eq . 
The DM clumps of mass M can be formed from density 
fluctuations of a different peak-height v = S cq /cr cq (M), 
where a cq (M) is a fluctuation dispersion on a mass- 
scale M at the time i e q- A mean internal density of 
the clump p is fixed at the time of the clump formation 
and according to [34| is p — np cq [va cq (M )/8 c ] 3 , where 
S c = 3(12tt) 2 / 3 /20 ~ 1.686. 

A tidal destruction of clumps is a complicated process 
and depends on many factors: the formation history of 
clumps, host density profile, the existence of another sub- 
structures inside the host, orbital parameters of individ- 
ual clumps in the hosts, etc. Only in numerical simula- 
tions, all these factors can be taken into account properly. 
The first such simulation in the small-scale region was 
produced in [l9| ■ We use a simplified analytical approach 
by parameterizing the energy gains in tidal interactions 
by the number of tidal shocks per dynamical time in the 
hosts. Using the model [35j for tidal heating, we deter- 
mine the survival time T, i.e. time of tidal destruction, 
for a chosen small-scale clump due to the tidal heating 
inside of a host clump with larger mass. During the dy- 
namical time idyn — 0.5(Gp/ l ) -1 / 2 , where pn is a mean 
internal density of the host, the chosen small-scale clump 
may belong to several successively destructed hosts. A 
clump trajectory in the host experiences successive turns 
accompanied by the "tidal shocks" [35[. Similar shocks 
come from interactions with other substructure, and in 
general due to any varying gravitational field. For the 
considered small-scale clump with a mass M and radius 
R, the corresponding internal energy increase after a sin- 
gle tidal shock is 

47T 

AE~ — nGphMR 2 , (38) 

where a numerical factor 71 ~ 1. Let us denote by 72 
the number of tidal shocks per dynamical time idyn- The 
corresponding rate of internal energy growth for a clump 
is E = 72AEytdyn- A clump is destroyed in the host if 



its internal energy increase due to tidal shocks exceeds a 
total energy \E\ ~ GM 2 /2R. As a result, for a typical 
time T = T(p, pt) of the tidal destruction of a small-scale 
clump with density p inside a more massive host with a 
density pn we obtain 

T-\p,p h ) = ||| * 4 7l72 G 1 / 2 ^ /2 p- 1 . (39) 

It turns out that a resulting mass function of small-scale 
clumps (see this Sec. below) depends rather weakly on 
the value of the product 7172- 

During its lifetime, a small-scale clump can stay in 
many host clumps of larger mass. After tidal disruption 
of the first lightest host, a small-scale clump becomes 
a constituent part of a larger one, etc. The process of 
hierarchical transition of a small-scale clump from one 
host to another occurs almost continuously in time up 
to the final host formation, where the tidal interaction 
becomes inefficient. The probability of clump survival, 
determined as a fraction of the clumps with mass M sur- 
viving the tidal destruction in hierarchical clustering, is 
given by the exponential function e~ J with 

j^y^L_ (40) 

Here, Ath is a difference of formation times th for two 
successive hosts, and summation goes over all clumps 
of intermediate mass-scales, which successively host the 
considered small-scale clump of a mass M. Changing the 
summation by integration in (|40p we obtain 

*/ 

7I s f dt h p x - p f p x t 2 

J{P,Pf) = 7f7 7-7 ~ 7 — ~7T2> 

J T{p,p h ) p p tf 

ti 

(41) 

where 

7 = 2 7l7 2 K 1 / 2 G 1 / 2 p^ oq ~ 14( 7l72 /3), (42) 

and t, ti, tf, p, pi and pf are, respectively, the formation 
times and internal densities of the considered clump and 
of its first and final hosts. One may see from Eq. (|4"Tj) 
that the first host provides a major contribution to the 
tidal destruction of the considered small-scale clump, es- 
pecially if the first host density pi is close to p, and con- 
sequently e~ J -C 1. Therefore, Eq. (|4Tj) gives a qualita- 
tively correct description of the tidal destruction. How- 
ever, in the more detailed approach the dependence of 7 
on another parameters is possible to take into account. 
As reasonable estimate, we will use the ansatz given by 
Eq. (|4"Tj) for further calculation of mass function. 

Now we need to track the number of clumps M (orig- 
inated from the density peak u) which enter some larger 
host during time intervals Atj around each ti beginning 
from the time t of clump formation. A mass function 
of small-scale clumps (i.e., a differential mass fraction of 
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DM in the form of clumps survived in hierarchical clus- 
tering) can be expressed as 



dM 



dM dv ■ 



(43) 



dti 



dM dh 



-J(t,ti) 



In this expression, to is the Universe age and F{M, t) is a 
mass fraction of unconfined clumps (i.e., clumps, not be- 
longing to more massive hosts) with a mass smaller than 
M at time t. According to [3J], the mass fraction of un- 
confined clumps is F(M,t) = erf (S c /[V2a cq (M)D(t)}) , 
where erf (x) is the error-function and D(t) is the growth 
factor normalized by D(t Qq ) = 1. An upper limit of inte- 
gration to in Eq. (|43[) is not crucial and may be extrapo- 
lated to infinity because a main contribution to the tidal 
destruction of clumps is provided by the early formed 
hosts at the beginning of the hierarchical clustering. 

Two processes are responsible for time evolution of the 
fraction d 2 F / (dMdt) for unconfined clumps in the mass 
interval dM: (i) the formation of new clumps and (ii) 
the capture of smaller clumps into the larger ones. Both 
these processes are equally efficient at the time when 
d 2 F/(8Mdt) — 0. To take into account the confined 
clumps (i.e., clumps in the hosts) we need only the second 
process (ii) for the fraction dF(M,t)/dM. Nevertheless, 
in Eq. (|l3"|l. which it is used the fraction 8F(M,t)/dM 
which depends on both processes. This is not accurate at 
a typical formation time of a clump with a mass M, when 
clump density is comparable with the density of hosts. 
Fortunately, for this time the exponent in Eq. (|43[) is very 
small, e _J <§; 1, as it can be seen from (|4"Tj) and AS)) . Re- 
spectively, an uncertain contribution from the process (i) 
to the integral (j4"3"|) is also very small. Meanwhile, only 
process (ii) dominates in the integration region where 
the exponent e~ J is not small. For this reason, Eq. (j4"3"|) 
provides a suitable approximation for the mass fraction 
of clumps survived in the hierarchical clustering. The 
characteristic epoch t* of the clumps M formation can 
be estimated from the equation o- eq (M)D(t*) ~ 5 C . If 
one considers, the times t ^> t*, then the exponents 
exp{— (5 2 /[2cr 2 q Z) 2 (£)]} can be putted approximately to 
unity for simplification of integration in (|43[) . 

Finally, we transform the distribution function (|4"3"|) to 
the following form: 



M 



^ e - 2/2 /i(7) dl ° S ^ (M) dM, 
\/2n 



dM 



(44) 
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FIG. 3: The function from (05 



/i(7) varies rather slowly in the interesting interval of 
14 < 7 < 40, and one may use fi(j) — 0.2 — 0.3. 

Physically, the first factor v in (|4"4")l corresponds to a 
more effective survival of high-density clumps (i.e., with 
large values of v) with respect to the low-density ones 
(with small values of v). Integrating Eq. (|4~4")) over v : we 
obtain 



, dM , dM 
6nt ^ 0.02 n + 3 . 

^ M V ' M 



(46) 



An effective power-law index n in Eq. (|46[) is determined 
as n = —3(1 + 2<91ogcr cq (M)/<91og M) and depends very 
weakly on M, Equation ([46]) implies that for suitable 
values of n only a small fraction of clumps, about 0.1 — 
0.5 %, survives the stage of hierarchical tidal destruction 
in the each logarithmic mass interval AlogAf ~ 1. It 
must be stressed that a physical meaning of the survived 
clump distribution function £ dM /M is different from the 
similar one for the unconfined clumps, given by the Press- 
Schechter mass function dF/dM. 

The simple M _1 shape of the mass function (|4"6"|) is in 
very good a gree ment with the corresponding numerical 
simulations [19(, but our normalization factor is a few 
times smaller. One also can see a reasonable agreement 
between the extrapolation of our calculations and the 
corresponding numerical simulations of the large-scale 
clumps with M > 10 6 M Q (for a comparison see Q). 
The obtained mass function (|46[) is further transformed 
in the process of tidal destructions of clumps by stars in 
the Galaxy (see previous sections). 



VI. MODIFIED DISTRIBUTION FUNCTION 
OF CLUMPS 



where 



A(7) 



2[T(l/3)-r(l/3,7)] 
3v / 2^7 1/3 



(45) 



r(l/3) and r(l/3, 7) are the Euler gamma-function and 
incomplete gamma-function, respectively. The function 
131) is shown in Fig. [31 It is seen in this figure that 



In this section we calculate the modified mass function 
for the small-scale clumps in the Galaxy taking into ac- 
count clump mass loss instead of the clump destruction 
considered in 0, Q3 • 

According to theoretical model [i[ and numerical sim- 
ulations [l9|], a differential number density of small- 
scale clumps in the comoving frame in the Universe is 
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FIG. 4: Numerically calculated modified mass function of 
clump remnants for galactocentric distances 3 and 8.5 kpc. 
The solid curve shows the initial mass function. 



n(M) dM oc dM/M 2 . This distribution is shown in Fig.[4] 
by the solid line. The damping of small-scale pertur- 
bations with M < M m i n provides an additional factor 
exp[— (M/M m i n ) 2 / 3 ] responsible for the fading of distri- 
bution at small M. The result of the numerical simu- 
lations can be expressed in the form of a differen- 
tial mass fraction of the DM clumps in the Galactic halo 
f(M) dM ~ K(dM/M), where k~8.3x 1CT 3 . The ana- 
lytical estimation (|4"rj|) gives approximately k~4x 10~ 3 
for the mass interval lCr 6 Af Q < M < IM Q . The dis- 
crepancy by the factor ~ 2 may be attributed to the ap- 
proximate nature of our approach as well as to the well 
known additional factor 2 in the original Press-Schechter 
derivation of the mass function. In the later case one 
must simply multiply equation (|43[) by factor 2. To clar- 
ify this discrepancy the more sophisticated calculations 
are necessary. 

As it was described earlier, we consider DM clumps 
originated from 2a density peaks. Therefore, in our ap- 
proximation the density of clumps and their distribution 
depends only on one parameter M . In general, the distri- 
bution of DM clumps depends on the pair of parameters, 
e.g. mass and radius, mass and velocity dispersion or 
mass and peak- high v as in the distribution (|44]) . Mean- 
while the authors of numerical simulations do not present 
a general distribution of clumps over two parameters. 
The general distribution of clumps can be in principle 
extracted from simulations and is very requested for fur- 
ther investigations of DM dumpiness. 

By using the formalism of Sec. El we derive the mass 
distribution of the clump remnants in dependence of the 
initial masses Mj of clumps. To do this, we calculate nu- 
merically the value of the mass fi of the clump remnant in 
dependence of the initial mass Mi for separate elements 
ApA^Ae in the parameter space in ([5rj]) . Then for fixed 
intervals A/j, of values of /i, we provide the summation 
of the weights of distribution function, which is given by 
Eq. (|36|) without symbols of integration and 0-function. 
By using the derived ^-distribution, we transform the ini- 
tial (cosmological) mass function of clumps to the final 



(nowadays) mass function in the halo at the present mo- 
ment. This final mass function is shown in Fig. 0] for two 
distances from the Galactic center. We supposed in nu- 
merical calculations that a core radius is very small and 
all masses of remnants are admissible. With a finite core 
size, the final mass function has a cutoff near the cores 
mass of clump with a minimal mass M m j n . The adiabatic 
correction leads to the accumulation of remnants of some 
mass corresponding to violation of momentum approxi- 
mation. One can see from Fig. U that clump remnants 
exist below the M m i n . Deep in the bulge (very near to the 
Galactic center) the clump remnants are more numerous 
because of intensive destructions of clumps in the dense 
stellar environment in comparison with the rarefied one 
in the halo. The main contribution to the low-mass tail 
of the mass function of remnants comes from the clumps 
with the near-disk orbits where the destructions are more 
efficient. 

The another important point is an efficient destruction 
of clumps with orbits confined inside the stellar bulge. 
Nevertheless, a number density of clumps inside the bulge 
is nonzero because a major part of clumps have orbits 
extending far beyond the bulge, These "transit" clumps 
spend only a small part of their orbital time traversing 
the bulge and survive the tidal destruction. 

VII. AMPLIFICATION OF ANNIHILATION 
SIGNAL 

A local annihilation rate is proportional to the square 
of the DM particle number density. A number density of 
DM particles in clump is much large than a correspond- 
ing number density of the diffuse (not clumped) compo- 
nent of DM. For this reason, an annihilation signal from 
even a small fraction of DM clumps can dominate over 
an annihilation signal from the diffuse component of DM 
in the halo. In this section, we calculate the amplifi- 
cation (or "boosting") of an annihilation signal due to 
the presence of the survived DM clump remnants in the 
Galactic halo. We consider here the Harrison- Zeldovich 
initial perturbation spectrum with power index n p = 1 as 
a representative example. The value of n p is not exactly 
fixed by the current observations of cosmic microwave 
background (CMB) anisotropy. In the case of n p < 1, 
the DM clumps are less dense, and a corresponding am- 
plification of annihilation signal would be rather small 

The gamma-ray flux from the annihilation of the dif- 
fuse distribution (fT5)) of DM in the halo is proportional 
to 

w(C) 

l n = J P 2 H (Odx, (47) 
o 

where the integration is over r goes along the line of sight, 
£(C r ) = ( r2 + r © — 2rr cosC) 1 / 2 is the distance to the 
Galactic center, r max (£) = (R^ — r| sin 2 C) 1 ^ 2 + r o cos£ 
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FIG. 5: The annihilation signal (|48[) (upper curve) as a func- 
tion of the angle £ between the line of observation and the 
direction to the Galactic center. For comparison the annihi- 
lation signal is also shown (by the bottom curve) from the 
Galactic halo without DM clumps (|47[) . The values of both 
integrals (|48|l and ()47[) are multiplied by a factor of 10 48 . 

is a distance to the external halo border, £ is an angle 
between the line of observation and the direction to the 
Galactic center, i?H is a virial radius of the Galactic halo, 
r Q = 8.5 kpc is the distance between the Sun and the 
Galactic center. The corresponding signal from annihi- 
lations of DM in clumps is proportional to the quantity 

a 

rw(C) 

I cl = S J dx J f(M)dMpp H (£)P((,p), (48) 

M min 

where p(M) is the mean density of the clump. The func- 
tion S depends on the clump density profile and core 
radius of the clump 4], and we use S ~ 14.5 as a rep- 
resentative example. The observed amplification of the 
annihilation signal is defined as r](Q = (I c i + In) /In is 
shown in Fig. [6] for the case x c = 0.1. It tends to unity 
at £ — * because of the divergent form of the halo pro- 
file (fT5")l . The annihilation of diffuse DM prevails over 
signal from clumps at the the Galactic center. The r](£) 
very slightly depends on x c , and corresponding graphs 
for x c < 0.1 are almost indistinguishable from the one in 
Fig. [5l This is because the observed signal is obtained 
by integration along the line of sight and the effect of the 
clump's destruction at the Galactic center is masked by 
the signal from another regions of the halo. 

This amplification of an annihilation signal is often 
called a "boost-factor" . A boost-factor of the order of 
10 is required for interpretation of the observed EGRET 
gamma-ray excess as a possible signature of DM neu- 
tralino annihilation 1361. 



VIII. CONCLUSION 

In Q it was found that almost all small-scale clumps 
in the Galaxy are destructed by tidal interactions with 



FIG. 6: The amplification of the annihilation signal (J c i + 
/h)//h as function of the angle between the line of observation 
and the direction to the Galactic center, where fluxes are given 
by flU} and (T351) . 



stars and transformed into "ministreams" of DM. The 
properties of these ministreams may be important for 
the direct detection of DM particles because DM parti- 
cles in streams arrive anisotropically from several discrete 
directions. In this work, we demonstrate that the cores of 
clumps (or clump remnants) survive in general during the 
tidal destruction by stars in the Galaxy. Although their 
outer shells are stripped and produce the ministreams 
of DM, the central cores are protected by the adiabatic 
invariant and survived as the sources of annihilation sig- 
nals. This conclusion depends crucially on the unknown 
sizes of the cores: the smaller cores are more protected 
because DM particles there have higher orbital frequen- 
cies and therefore the larger the adiabatic parameter. 

Despite the small survival probability of clumps dur- 
ing early stage of hierarchial clustering, they provide the 
major contribution to the annihilation signal (in compar- 
ison with the unclumpy DM). The amplification (boost- 
factor) can reach 10 2 or even 10 3 depending on the ini- 
tial perturbation spectrum and minimum mass of clumps. 
This boost-factor must be included in calculations of the 
annihilation signals. Some promising interpretations of 
observations and calculations of annihilation signal from 
the Galactic halo require this boost-factor (see, e.g., [36|]). 
The discussed dense remnants survive the tidal destruc- 
tion and provide the enhancement of DM annihilation 
in the Galaxy. These remnants of DM clumps form the 
low-mass tail in the standard mass distribution of small- 
scale clumps extended much below M m i n of the standard 
distribution. It does not mean of course the increasing of 
annihilation signal in comparison with the case without 
clump destruction. It only indicates that galactic clump 
destruction does not diminishes strongly the annihilation 
signal. 

The principle simplifying assumption of this work is 
that only the outer layers of clumps are subjected by the 
tidal stripping. The main difficulties in considering the 
full problem with the mass loss from inner layers are in 
the complicated dynamical reconstruction of clumps just 
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after tidal shocks. We believe that our approach provides 
a rather good result by the two reasons. First, the in- 
fluence of tidal forces depends on the system size, and, 
therefore, the outer layers are greatly subjected to tidal 
forces. Second, the adiabatic protection is more efficient 
in the inner part of the clump because of higher orbital 
frequencies here. In reality, we expect some expansion of 
clump and diminishing of its central density due to en- 
ergy deposit from tidal forces. It is very interesting task 
to clarify this process in future works. 

The numerical estimate of the boost-factor for DM par- 
ticle annihilation inside clumps is very model-dependent. 
It depends on nature of DM particles and on their in- 
teraction with ambient plasma. The important physical 
parameters, which affects the annihilation rate in clumps, 
are decoupling temperature Td and minimal mass M m i n 
in the clump mass distribution. The boost-factor in- 
creases strongly for small M m ; n . The minimal mass in 
standard calculations is determined by the escape of DM 
particles from a growing fluctuation due to, e.g., diffu- 
sion, free streaming or Silk effect. Uncertainties in the 
calculated values of Td and M m i n are discussed in Ap- 
pendix IA II For the lightest neutralino as a DM par- 
ticle, assuming it to be the pure bino B, one can see 
from Table III a huge difference in M m ; n caused by 
the variation of supersymmetry (SUSY) parameters m x 
and m. For these parameters, we use cosmologically al- 
lowed values from the benchmark scenarios of the work 
[37l ]. Moreover, inclusion of other neutralino composi- 
tions, e.g., mixed bino-Higgsino, the other allowed bench- 
mark scenarios with co-annihilation and focus-point re- 
gions, and some other modifications, may very consid- 
erably increase the allowed region of M m ; n values up to 
(3 x 10" 12 - 7 x lCr 4 )Mo ||]. Inclusion of the other 
particle candidates extends further this region. 

Another parameter variation which affects strongly the 
boost-factor is the spectral index of density perturba- 
tion rip (see [H). We conclude thus that the annihila- 
tion boost-factor (enhancement) even for neutralino has 
large uncertainties due to the difference in SUSY pa- 
pameters and spectral index n p . It can reach the factor 
10 4 and even more, The largest values of boost-factor can 
be already excluded by observations of indirect signal, 
since mSUGRA parameters can be fixed for this largest 
value. On the contrary, a tidal destruction of clumps 
in the Galaxy affects the annihilation boost-factor much 
weaker. 
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APPENDIX A: MINIMUM CLUMP MASS 

1. Uncertainties in minimal clump mass and 
decoupling temperature 

The low-mass cutoff of the clump mass-spectrum ac- 
companies the process of decoupling. It starts when DM 
particles coupled strongly with surrounding plasma in 
the growing density fluctuations. The smearing of the 
small-scale fluctuations is due to the collision damping 
occurring just before decoupling, in analogy with the Silk 
damping [39| ■ It occurs due to the diffusion of DM parti- 
cles from a growing fluctuation, and only the small-scale 
fluctuations can be destroyed by this process. The corre- 
sponding diffusive cutoff is very small. As coupling 
becomes weaker, the larger fluctuations are destroyed 
and M m ; n increases. One may expect that the largest 
value of M m j n is related to a free-streaming regime. How- 
ever, as recent calculations show [l(|, the largest M m i n 
is related to some friction between DM particles and 
cosmic plasma similar to the Silk damping. The pre- 
dicted minimal clump masses range from very low values, 
^Umin ~ 10~ 12 M© [40], produced by diffusive escape of 
DM particles, up to M m ; n ~ 10 _4 M Q , caused by acous- 
tics oscillations [42| and quasi-free-streaming with limited 
friction [icj |. 

The calculations of minimal clump mass Af m ; n and de- 
coupling temperature Td are determined by elastic scat- 
tering of DM particles off leptons I = [yhi^hi&R) in 
cosmic plasma. The uncertainties in cross-section very 
strongly influence the resulting values of M m ; n and 
In all works cited above, the lightest neutralino (x) in 
the form of pure bino (B) is assumed as a DM parti- 
cle, and x^-scattering occurs due to exchange by sleptons 
I = {vL,eL,e.R). 

The elastic cross-sections for and IrX scattering 
have been calculated in Q as 

fda\ a 2 m ^ 2 (1 + cos 6> cm ) (A] 

\dnj lLX 8cos 4 w (m 2 -™ 2 ) 2 > 

and 

(10, - >•(£), ■ if *>-*■* (A2) 

v / IrX v / IlX 

where w > m; is a cm. -energy of I, 6> cm is a scattering 
angle of I in c.m.-system, m x is a neutralino (bino) mass, 
fhjj and rhu are, respectively a mass of the left and right 
sfermions, and $w is the Weinberg angle. 

The values of T d and M min as cited in 0, 1, [M SI 
differ very much from each other, but a very big contri- 
bution to this difference comes from the differences in the 
used values for m x and m. To sec the difference, which 
must be attributed to the different damping mechanisms 
used in these works, we recalculated Td and M m - ln with 
the same values of m x and m, for which we used 100 
and 200 GeV respectively. The results are presented in 
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Reference: [41] 1} [4] 1} [6] 2) [42] 3) [10] 4) rh Td M min 

T d , MeV 28 26 25 20 22.6 100 GeV 200 GeV 22.6 MeV 8.4 x 1O" 6 M 

M min /M s 2.5 x 10" 7 1.7 x 10 -7 1.5 x 10 -6 1.3 x 10 -5 8.4 x 10 -6 100 GeV 1500 GeV 196 MeV 1.3 x 1O _8 M 

800 GeV 1600 GeV 305 MeV 3.5 x 1O- 9 M 

TABLE I: The values of decoupling temperature Td and min- 
imal clump mass Af min with m x = 100 GeV and rh = TABLE III: Values of T d and A/ min for the Bertschinger ^ 
200 GeV for different damping mechanisms: Tree-streaming, damping scenario and three benchmark scenarios [37| which 
'collision damping, 'acoustic oscillations, 'quasi-free- close to scenarios B ', E' and M' shown in the Table II. 
streaming with friction. 



scenario % &l £r v e , £y v T 

B' 95 188 117 167 167 

E' 112 1543 1534 1539 1532 

M' 794 1660 1312 1648 1492 

TABLE II: Selected benchmark scenarios from 37]. The 
masses of particles are given in GeV. 



Table I. We did not include there the work [4(| because 
a pure diffusive damping results in too low a value of 
-Wmin- From the Table I one can see a reasonable 
agreement in values of Td and M m j n , and a successive 
increasing of M m - m from 2.5 x 1O _7 M0 for free-streaming 
to ~ 1O -5 M0 for oscillation damping and quasi-free- 
streaming with friction. 



2. Uncertainties in SUSY parameters 

We shall consider now the range of predictions for dif- 
ferent values of SUSY parameters allowed in cosmology. 
For this aim, we shall use the SUSY benchmark scenar- 
ios from the work [371 ] . which agree with the Wilkinson 
Microwave Anisotropy Probe (WMAP) and other cos- 
mological data. These benchmark scenarios are obtained 
within mSUGRA model with universal parameters at the 
grand unified theory scale: mo (the universal scalar soft 
breaking mass), (the universal gaugino soft break- 

ing mass), A (the universal cubic soft breaking terms) 
and tan/3 (the ratio of two Higgs v.e.v.'s). The LEP and 
b — > S7 constraints are imposed. The resulting relic den- 
sity of neutralinos from these scenarios is in agreement 
with the WMAP data or can be obtained with small 
changes of m and mi/2- I n Table II we display three 
benchmark scenarios from [37|]- The scenario B' gives 
the lower value m x w 100 GeV and m close to 200 GeV, 
which we discussed above. Scenario M' gives the highest 
value m x w 800 GeV and m w 1600 GeV. Respectively, 
scenario E' gives the intermediate value m x ss 110 GeV 
and m « 1500 GeV, similar to those we used in 

To illustrate the uncertainties in Td and M m j n due to 
uncertainties in m x and rh (in the simplifying assumption 
that mc, = Trie, = mg H ) we choose the calculations of 
Bertschinger [lOj in quasi-free-streaming scenario with 
friction, which seem to be at present the most detailed 



ones. Wc use the Bertschinger formulaes 

M min = 7.59 x 10- 3 C 3 / 4 ( ^ Q X ^ ) ' M Q , (A4) 
with a dimensionless constant 

C = 256 (G F m 2 w ) 2 - l) £>i + 4), (A5) 

where Gf is the Fermi coupling constant, 6l and cj, are 
left and right chiral vertices, and tow, to, and m x are, re- 
spectively, the masses of the W boson Gprn^y — 0.0754, 
the slepton, the neutralino and the number of freedom at 
the decoupling epoch g» = 43/4. (Our own calculations 
in Q of C, which is related with the square of the matrix 
element for I + x — ► I + X scattering, differ from (|A4I) by 
a factor 1.6.) As a result we obtain for the benchmark 
scenarios which approximately coincide with model B' 
(minimum m x and rh), E' (minimum m x and large rh) 
and M' (very large m x and to) the values of Td and M m - ln 
listed in Table III. The predicted range of parameters for 
M min from this Table: 3.5 x 10" 9 -8.4x 1O" 6 M is not ro- 
bust at all. It is obtained within mSUGRA assumptions 
about possible universality of SUSY parameters too, mi/2 
and Aq. Lifting the universality restriction, the mass of 
the neutralino can increase up to the TeV range scale 
(though to v > 200 GeV needs a fine-tuning less than 1% 
in SUSY [43| or decreased down to a few GeV 0). 

In the numerical predictions above, we limited our- 
selves by rather restrictive assumptions on the mSUGRA 
model. The most important of them are assumption 
that the neutralino is a pure bino state and a choice 
of cosmologically allowed benchmark scenario. The de- 
tailed analysis made in [38[ showed that allowed parame- 
ters of mSUGRA result in much wider possibilities, e.g., 
neutralino as mixed bino-Higgsino and the other bench- 
mark scenarios. These possibilities are considered un- 
der WMAP cosmological constraints and a condition 
of producing the corresponding DM density for each 
set of mSUGRA parameters. The considered modifica- 
tions allow new channels of neutralino interactions with 
ordinary particles, e.g., the exchange by Z-boson, co- 
annihilation and resonances in neutralino-fermion scat- 
tering. It results in a wide range (many orders of mag- 
nitude) of scattering cross-sections, and, respectively, in 
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a wide range of decoupling temperature from 5 MeV to 
3 GeV. The corresponding range of M m [ n is given by 
(3 x 1CT 12 - 7 x 1CT 4 ) M Q . The authors consider also 
the Kaluza-Klein particle as a DM candidate. 

The small-scale mass of M m j n results in the large den- 
sity of DM clumps, and thus in a much stronger annihi- 
lation signal from the Galactic halo. For typical values 
of power-index of perturbation spectrum (from CMB ob- 



servations) the small-scale mass of M m i n results in the 
large density of DM clumps, and thus in a much stronger 
annihilation signal from the Galactic halo. However, the 
dependance of a mean clump density on the clump mass 
is rather weak due to the nearly flat form of the pertur- 
bation spectrum at small scales. The crucial factor for 
the amplification of the annihilation signal by clumps is 
the value of the perturbation power-law index n p . 
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